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' Thermodynamic quantities of Coulomb plasmas consisting of point-like ions immersed in a com- 

pressible, polarizable electron background are calculated for ion charges Z = 1 to 26 and for a 
wide domain of plasma parameters ranging from the Debye-Hiickel limit to the crystallization point 
and from the region of nondegenerate to fully degenerate nonrelativistic or relativistic electrons. 
The calculations are based on the linear-response theory for the electron-ion interaction, includ- 
' ing the local-field corrections in the electronic dielectric function. The thermodynamic quantities 

are calculated in the framework of the A'^-body hypernetted-chain equations and fitted by analytic 
expressions. We present also accurate analytic approximations for the free energy of the ideal elec- 
i-^ ' tron gas at arbitrary degeneracy and relativity and for the excess free energy of the one-component 

plasma of ions (OCP) derived from Monte Carlo simulations. The extension to multi-ionic mixtures 
is discussed within the framework of the linear mixing rule. These formulae provide a completely 
^ , analytic, accurate description of the thermodynamic quantities of fully ionized electron-ion Coulomb 

■ plasmas, a useful tool for various applications from liquid state theory to dense stellar matter. 

Ph; PACS numbers: 52.25.Kn, 05.70.Ce 

o 

I. INTRODUCTION 

Electron-ion plasmas (EIP) consisting of different species of point-like ions (charge Z^e, mass = Ai a.m.u.) and 
electrons (— e, rrie) are encountered in numerous physical and astrophysical situations like, e.g., inertially confined 
laboratory plasmas, liquid metals, stellar and planetary interiors, supernova explosions, etc. |l|]. Full ionization is 
^ ' reached either at high temperatures T and low densities p (thermal ionization) or at high enough densities p (pressure 
ionization). Even when these conditions are not satisfied, the approximation of full ionization is useful for calculations 
^~L^ in the mean-ion approximation, in which the mean ion charge corresponds to its partial ionization stage. On the 
J — ] other hand, the free energy of fully ionized EIP provides the reference system for models aimed at describing the 
• thermodynamic properties of partially ionized plasmas In this paper, we present a completely analytic model 
OO [ for the free energy of EIP, based on detailed numerical calculations for different ionic species Z over a wide range of 
0^ . density and temperature. We first focus on the two-component plasma (TCP), consisting of electrons and a single 

' species of ions. Extension to ionic mixtures is considered in Sect. VI. 
O , The Coulomb plasmas can be characterized by the electron coupling parameter Pg and the density parameter Ts, 



- ' — \ 

X 



Te = /Se^/fle, = ae/a_B, (1) 

^ ^' where (3 = {ksT)^^ is the inverse thermodynamic temperature, ks is the Boltzmann constant, = {^T:ne)~^^'^ 
^ ' measures the mean inter-electron distance, Ug is the electron number density, and as — /(mee^) is the Bohr radius. 
These parameters can be evaluated as Pe « (2.693 x 10^ K/T) 7124'^ = (23.2 eV/ZcsT) 7124^ and Ts ~ 1.172 7^24^^'^, where 
7124 = rie/lO^^ cm~'^ « (p/1.6605 g cm~^){Z) / {A) . Here and hereafter, (A) = "^-riiXi /J2i^i denotes the average 
over all ions and rii the number density of ions of zth species. 
The ion coupling parameter of the TCP is 

P, =/3(Ze)Va, = PeZ5/3, (2) 

where a; is the mean interionic distance {ai = a^Z^^^ due to the electroneutrality condition rig = UiZ). In multicom- 
ponent plasmas (MCP), it may be useful to define P, = Pe(^^/^). 
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The degeneracy parameter 9 and the relativity parameter x are defined respectively as: 

6 = T/Tf, X = pF/{mec), (3) 

where Tp — (meC^/kB) [\/T+~i? — 1] w (5.93 x 10^ K) [\/T+~e^— 1] is the Fermi temperature, c is the speed of light, 
and pf = 7i(37r^?ie)^/'^ is the zero-temperatm'e Fermi momentum of electrons. To estimate 9 and x, it is useful to note 
that 

_ /97ry/^ a _ 0.014 ^ f (Z) p 
^ ~ V 4 / rs V (A) 106 g cm-3 

Q^^^eTsT^ 0^ 0.543— ata;«l, (5) 

where a = 1/137.036 is the fine-structure constant. 

Various asymptotic expansions, interpolation formulae and large tables have been derived over the past for the 
thermodynamic functions of free fermions (see Refs. and references therein). In this paper, we present ana- 
lytic expressions for the thermodynamic quantities of free fermions for arbitrary degeneracy and relativity, 9 and x. 
Secondly, we propose simple and accurate analytic approximations for the nonideal internal and free energies of the 
classical one-component plasma (OCP), which take into account the most recent hypernetted-chain (HNC) and Monte 
Carlo (MC) calculations by DeWitt et al. Q (DWSC) in the strong-coupling regime. Third, we consider the electron 
screening effects on the thermodynamic properties of the TCP. We employ a computational HNC scheme based on 
the linear screening theory with local field corrections, taking into account the finite temperature (finite 9) effects. 
The numerical calculations have been performed over a wide range of Z, F^, and r^, and interpolated by a simple 
analytic formula, which recovers the Debye-Hiickel (DH) limit for the TCP at Ti ^ 1 and the Thomas-Fermi limit at 
large F^ and Z. 




II. SUMMARY OF THE MODEL 



Consider the Helmholtz free energy F, internal energy U, and pressure P of a TCP of iV^ ions and Ne electrons in 
the volume V. The total free energy Ftot can be written as the sum of three terms, 

Ftot = ^^/;^+^^/;^+^^cx, (6) 

(i e) . 

where F^^' denote the ideal free energy of ions and electrons, respectively, and -Fox is the excess free energy arising 
from interactions. 

In this paper, we restrict ourselves to conditions where the ions behave classically, which is the case in most 
astrophysical situations. Quantum corrections for ions which can be important in the ultradcnse matter of white 
dwarf interiors, neutron stars, and supernova cores have been considered, e.g., in Refs. [^[0]. Thus Fj^' is given by 

the Maxwell-Boltzmann expression. For Fj^^' , we use the well-known expressions of the thermodynamic functions of 
the perfect gas of fermions (which may be degenerate and relativistic) through the generalized Fermi-Dirac integrals. 

To calculate Fcx, we follow the model developed by Chabrier Q for fully ionized EIPs. As long as the ion-electron 
interaction is weak compared to the kinetic energy of the electrons, Ze'^/ue <C UbTf^ this interaction can be treated 
within the linear screening theory. Under these conditions, the exact Hamiltonian of the TCP can be separated out 
exactly into a Hamiltonian for the electron-screened ionic fluid and a Hamiltonian for a rigid electron background, 
the so-called "jellium" Hamiltonian 

H = H"^ + iJe (7) 

with 

where Ki is the ionic kinetic (translational) term, is the Fourier component of the ionic microscopic density and 
e{k) is the static screening function of the electron fluid to be discussed below. The Hamiltonian H'^^ characterizes 
the electron-screened ion fluid with the interparticle potential whose Fourier transform is 



e(k) 



N,, 



(8) 



2 



(9) 



which is the sum of the bare ionic potential and the induced polarization potential. 

The ion-ion (ii) and the ion-electron (ie) Coulomb interactions can thus be separated from the exchange-correlation 
contribution in the electron fluid (ee). The excess part of the free energy can then be written as i^ex = Fee+Fu+Fie] 
the quantities labeled ie will be referred to as electron- screening quantities. It is convenient to consider dimensionless 
quantities fee = PFee/Ne and fa^ie = PFn^ie/Ni. Then 

where Xi^e = Ni^e/N denote the number fraction of ions and electrons, respectively, and N = Ni + Ne is the total 
number of particles. In the same way we define Uee = fiUee/Ne and uu^ie = pUu^ie/Ni. The excess free energy can be 
obtained from the internal energy by integration: 

/cx(r,r.) = ^'^^;i^dr. (11) 

For /ee, we have adopted the interpolation formula of Ichimaru ct al. [Tl| ] (hereafter IIT), consistent with numerical 
results obtained by different authors. 

For /ii, which corresponds to the well-known OCP model, that implies the rigid electron background (e(fc) = 1), 
we present an analytic interpolation between the MC results [|| at Ti > 1 and the DH limit and Abe correction at 

r. <o.i. 

The ion-electron interactions are calculated numerically as in Ref. In this approach, the bare Coulomb potential 
in the expression for the electrostatic energy is replaced by the potential statistically screened by the electrons (|^), 
and the HNC approximation is used to calculate the thermodynamic functions of the system. This model, originally 
applied to nonrelativistic hydrogen plasmas, is now extended to the case of arbitrary Z and x. In the nonrelativistic 
case {x 1), the dielectric function e{k) is the finite-temperature Lindhard function modified with the local field 
correction arising from electron correlation effects. 

At very high density, x ^ 1, the electrons become relativistic. At such densities, the electron correlation effects are 
completely negligible. The finite-temperature effects {9 ^ 0) may give an appreciable contribution to the screening 
part of the free energy jie only at extremely high temperatures, where the nonideality of the gas has no significance. 
Thus we use the Jancovici |12(| zero-temperature dielectric function in the relativistic regime. 

The correlation functions and thermodynamic quantities for the electron-screened ionic fluid are obtained within 
the framework of the HNC equations. The validity of the HNC theory for the Coulomb systems has been assessed 
by several authors by comparison with lengthy MC simulations. The HNC approximation consists of neglecting the 
contribution of the so-called bridge diagrams, which involves an infinite series of multiple integrals, in the iV-body 
general diagrammatic resummations p^ . The long-range part of the direct correlation function c(r) calculated within 
the HNC approximation is exactly canceled by —V{r)/kT, so that the pair correlation function g{r) is of much shorter 
range than the Coulomb potential V{r) |T^. This is a required condition for Coulomb systems, because of the perfect 
screening condition. This property of the HNC theory makes it particularly suitable for such long-range systems. The 
differences on the free energy, the internal energy and the pressure are at most of the order of 1% (see, e.g., Refs. 
The difference is due to the lack of bridge functions in the HNC theory. 



III. IDEAL PART OF THE FREE ENERGY 

The ideal free energy of nonrelativistic classical ions, neglecting their spin statistics, reads p4[ : 

i^/;^ = N^ksT [Hn^XD - 1] , (12) 
where Xi = (27r/3?i^/mi)^/^ is the thermal wavelength of ions. For electrons, we use the identity 

F^:^=Ne^^li'-pl:^v. (13) 

Here, /i^^-* is the chemical potential (in which we do not include the rest energy meC^) and P^-^'^ is the pressure of the 
ideal Fermi gas. The pressure and number density, in turn, are functions of fi and T: 
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^ ' (4/2(X,T) + -/5/2(x,r)J, 

(-^1/2 (X,t) +T"/3/2(X,t)) 



id 



where r = (/SmeC^)-! = T/5.93 x 10^ K, x = /3Mid\ 



37r2;i3^5/2 

^/2(me//3)3/2 



and 



x'' ^1 + Ta;/2 
exp(a; — x) + 1 



da: 



(14) 
(15) 

(16) 



is the generalized Fermi-Dirac integral. 

In the hmit r — > 0, the Fermi-Dirac integrals reduce to the usual nonrelativistic Fermi integrals Iu{x)i which can be 
calculated using the highly accurate Fade approximations presented by Antia |T5| . The chemical potential is obtained 
from the relationship 

X = Xi/2(20-3/V3), (17) 

where Xy is the inverse Fermi integral, also fitted with high accuracy by Antia p^ . 

The accuracy of the nonrelativistic formulae decreases rapidly at T > 10^ K. Blinnikov et al. |^ have presented a 
number of approximations and asymptotic expansions of the relativistic thermodynamic functions of the ideal electron 
gas. We have selected those of their fitting formulae that are most accurate at low and moderate x ^-iid supplemented 
them with asymptotic expansions at high x to obtain an approximation which is accurate at any for each of the 
Fermi integrals Iy{Xi t) with v 



2 ' 2 ' 2 



4+1/2 (X,t) =^C, 



l + xfV/2 



1=1 

5 

E 

i=l 



exp(-xf^) + exp(-x) 

hiX. 



(X < 0.6), 



t4^^?^^4^^ + «.te+x)'=+^/Vl + te+x)r/2 
exp(xa;j - x) 



1 

(0.6 < X < 14), 

^2 



Fkix,r) + YX 



TT' j^A: + l/2+(fc + l)xT/2 



R 



(X > 14), 



where R = Vx(l + X^/2), 



^o(x, r)^ix + r-')R/2 ~ (2r)-3/2 in(l + rx + V2^ R), 
Fi{x,T)^ {2Ry3~Fo{x,r)) /t, 
F2{x,t) ^ {2xR^ ^5F,{x,t)) /{At). 



(18) 



(19) 
(20) 



(21) 
(22) 
(23) 



If XT ^ Ij the functions Fk{x,T) should be replaced by their nonrelativistic limits, x^^^^"^ + 3/2). The constants 
, Xi, ^i, hi, and Vi are adopted from Ref. |3| and listed in Tabic ^ The relative error of the approximation 



([r8|)-(p0|) does not exceed 0.2% at r < 10^ (any x), being typically a few parts in 10**. 

TABLE I. Parameters of Eqs. (|l|) and (^. The powers of 10 are given in square brackets 



i 


1 


2 


3 




4 




5 






0.37045057 


0.41258437 


9.777982 [- 


2] 


5.3734153 


-3] 


3.8746281 


-5] 




0.39603109 


0.69468795 


0.22322760 




1.5262934 


-2] 


1.3081939 


-4] 


xf> 


0.76934619 


1.7891437 


0.70754974 




5.6755672 


-2] 


5.5571480 


-4] 


0.43139881 


1.7597537 


4.1044654 




7.7467038 




13.457678 




xf' 


0.81763176 


2.4723339 


5.1160061 




9.0441465 




15.049882 




(2) 
Xi 

Xi 

hi 

Vi 


1.2558461 
7.265351 [-2] 
0.26356032 
3.818735 [-2] 
0.29505869 


3.2070406 
0.2694608 
1.4134031 
0.1256732 
0.32064856 


6.1239082 
0.533122 
3.5964258 
0.1986308 
7.3915570 [ 


-2] 


10.316126 
0.7868801 
7.0858100 
0.1976334 
3.6087389 


-3] 


16.597079 
0.9569313 
12.640801 
0.1065420 
2.3369894 


-5] 
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The chemical potential fi\^ can be obtained numerically from Eq. (p^), using Eqs. (p^)-(]20|). We have constructed 
also an analytic fit to x- 

Here T^n°'"'°i is given by the nonrelativistic formula (p^, and the coefficients qi are functions of 9: 

|(o'-l)-, 

12 + 80-^/^, 

2 e-'' + f.612e^ 

^ " 6.192 6I0 0944 c-9 ^ 5 g/o.egs ' 

The relative error Sx/x becomes infinite at % = 0. However, since thermodynamic quantities are expressed 
through X by virtue of thermal averaging of type of Eq. (|l^), a natural measure of the error is Sx/ max(|x|, 1) = 
(5/i/ max(|/i|, fc^T). The error thus lies within 0.4% for r > 1 and is smaller than 0.2% if r < 1 (any 9). Another 
measure of the accuracy is the relative difference between the densities calculated with the exact and fitted values 
of /i. This difference lies within 0.4% for r > 1 and within 0.f% for r < 1. 

This accuracy may not be sufficient for calculation of temperature derivatives of the electron-gas EOS (heat capacity, 
temperature exponent, etc.) in the regime of strong degeneracy (y 3> 1). In this case, however, Sommerfeld asymptotic 
expansions for these quantities may be used (see, e.g., Ref. fl^]). 

In this paper we do not consider the positrons, which are efficiently created at r > 1 (see Ref. |^ for description of 
the equilibrium electron-positron plasma). 



1 + t/{29) 



qi 



qiq^T 



1 + q2T 



(24) 



91 = 

92 = 



IV. OCP LIQUID OF CLASSICAL IONS 



Liquid and solid phases of the OCP have been studied extensively by various numerical methods, MC simulations 
or A^-body semianalytic theories like the HNC theory (see Refs. 11 1 for detailed reviews). All the thermodynamic 
functions of the OCP of classical ions in a uniform (rigid) electron background can be expressed as functions of the only 
parameter F^. The melting point of the OCP corresponds to P^ « 172, above which it forms a Coulomb crystal p8| . 
The most accurate MC and HNC results for the internal and free energies of the liquid OCP for 1 < Pi < 160 have 
been obtained recently by DWSC ||^ (see references therein for earlier results). The high precision of the calculations 
allowed the authors to investigate the tiny effects of nonadditivity of the excess energy of binary ionic mixtures, as 
will be commented in Sect. VI. 

DWSC have also derived a highly accurate analytic fit to the MC simulations of the internal energy of the OCP in 
the aforementioned P-range: 



= aV, + bVl 



c, 



(25) 



with a = -0.899126, h = 0.60712, c = -0.27998, and s ^ 0.321308. The maximum relative difference between 
calculated and fitted values reaches 17 parts in 10^ at P^ — 3.1748. 

Equation (|2^), however, does not apply to the weak-coupling region P^ < 1. At very small P^, the internal energy of 

the OCP must recover the well-known DH expression uu = — (\/3/2)P^^^ whereas at moderately small P^ this limit 
must include the Abe correction M: 



ln(3P,) 



7 1 



(26) 



where 7 — 0.57721 ... is the Euler's constant. 

We represent the internal energy of the ionic fluid (P^ ^ 170) by a simplified version of the fitting formula proposed 
by Hansen MM, 



3/2 



A, 



(27) 



where Ai and A2 are fitting parameters, and A3 ~ — V3/2 — A\l\fA^. We have found that the minimum relative 
difference between Eq. ( p7| ) and the MC results of DWSC |^ — smaller than 6 parts in 10'* — is obtained with 
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Ai = —0.9052 and A2 = 0.6322. This accuracy is sufficient for our present study since it is much better than the 
available numerical accuracy of the complementary contribution to the internal energy, Uig. As mentioned in Sect. ||, 
the HNC calculations of the sum ua + Uie ensure an accuracy of the order of 1%. 




-2 -1 1 2 



FIG. 1. Comparison of the fit (solid line) given by Eq. ( p7[ ) for the OCP internal energy with the DH and Abe asymptotic 
expansions at small Ti (dot-dash Unes) and with the DWSC results |^ at 1 < Fi < 160 (dots and circles). The dashed curve 
presents the interpolation of SB . 

Figure |l] presents a comparison of our interpolation formula ( pTf ) with the DH-Abe formulae, the MC results and 
fit (25) of DWSC 1^, and the interpolation proposed by Stolzmann and Blocker (hereafter SB) following Ebcling 
psf . Unlike SB, our Eq. ( p7| ) accurately reproduces Eq. ( p6| ) in the range F; ~ 0.01 — 0.1 and provides a smoother 
transition between the strong- (F^ > 1) and weak- (F^ <C 1) coupling regimes. 

Using Eqs. (|^) and (0) we obtain the Helmholtz free energy (cf. Ref. pQ]): 

/,,(F,) = ^1 yT:{A2 + T,)- A2 in(v/r,M2 + Vi + r.Ma) 

-H2A3 VTi"- arctan (^VTi") ■ (28) 

At Ti > 1, this formula gives fa which differs from the HNC calculations and fit of DWSC ||] by no more than 0.8%. 
This difference approximately coincides with that between the MC and HNC results for uu, therefore it should be 
attributed to the lack of the bridge functions in the HNC approximation (see Sect. ||). On the other hand, Eq. (|2S 
recovers the DH-Abe free energy with an error smaller than 0.6% at F^ < 0.1. 



V. ELECTRON FLUID 



The exchange and correlation effects in electron fiuid were studied by many authors. For instance, Tanaka et al. 
[24 1 calculated the interaction energy of the electron fluid at finite temperature in the Singwi-Tosi-Land-Sjolander 
[25 1 approximation and presented a fitting formula that reproduces their results as well as various exact limits with 
disgressions less than 1% (in particular, their formula incorporates the parametrization of the exchange energy by 
Perrot and Dharma-wardana |^). We adopt a modification of this formula given by HT [ pl| . 

The exchange-correlation free energy, fee, is obtained by integration from Eq. (|l l|) . It is important to note that 
Tanaka et al. |2^ give a fit to the interaction energy of the electron fluid but not to the thermodynamic internal 
energy (the quantities differ at flnite 9). This enabled Tanaka et al. to obtain fee by integration of their fitting formula 
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over Fe at constant (the integration of the internal energy would have to be performed at constant Tg). Note also 
that the results of IIT are nonrelativistic. 

More recently, SB |2|] proposed other parametrizations of the exchange and correlation free energies. At moderate 
r^, comparison of the formulae given by SB and IIT reveals only small differences, which do not exceed the uncertainty 
in the various numerical results found in the literature 
in the relativistic case. On the other hand, the SB fit reaches the classical OCP limit at large and moderate Fg 
with disgressions up to 4.4%, while the parametrization of IIT is several times more accurate in this limit. We shall 
use the IIT's formula hereafter. 



27 1 . Unlike IIT, SB evaluated the exchange energy aX 6 <\ 



VI. ELECTRON SCREENING 



A. Numerical calculations 



In order to calculate the screening contribution, we have employed the model of Ref. outlined in Sect. y. The 
HNC equations were solved numerically for the effective screened interionic potential (H) to obtain fa + /^e, uu + Wie, 
and Pa + Pie, and for the bare Coulomb potential to obtain fa, uu, and Pa. The difference represents the screening 
(«e) contribution to the thermodynamic quantities. 

The previous numerical results |^ have been obtained for the hydrogen plasma {Z — 1). We extend the calculations 
to different values of Z and a larger set of . Figure |2| shows the effective potentials for Z = 6 and Z = 26 at 
several values of rg, compared with the bare Coulomb potential and with V^^ in the zero-temperature {0 = 0) RPA 
approximation (no local field correction) . One can see that the latter approximation works well at the small value of 

= 0.0256 (lower panels), while it breaks down completely at = 1 (upper panels). 




FIG. 2. Effective ion- ion potentials at various approximations for Z — Q (left panels) and Z = 2Q (right), for two 
densities, = 1 (upper panels) and 0.0256 (lower ones) and three values of F^ — 1, 10, and 40. Solid lines represent 
finite-temperature V'^ including the local field correction, and dashes show the zero-temperature RPA approximation (the 
dashed and solid lines practically coincide on the lower left panel). Bare Coulomb potential is drawn by dots for comparison. 



The bulk of the calculations has been performed in the nonrelativistic approximation, for 13 ion charges from Z = 1 
to Z = 26 listed in the first column of Table |l[ at 10 values of the density parameter ranging from = 0.0256 
to Vg « 2, and, at each Z and r^, for several tens of values of the coupling parameter F^ which range from the DH 
limit at Fj = 0.001 to Fj ~ 200. As an example, calculated values of the normalized screening part of the free energy 
fie at Z = 6 are shown by filled circles in Fig. ^. Note that it is the account of the finite electron temperature in the 
dielectric function that allows us to reach the correct TCP DH limit at low values of F^ (see Ref. 1^). 
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FIG. 3. Nonrelativistic finite-temperature (filled circles) and relativistic zero-temperature (open triangles) calculated values 
of the screening part fie of the free energy of TCP for Z — 6, compared with the fit (j29|) (solid lines). The fit is also compared 
with the approximation of YS |l6| (dashed line) , which is valid at small and large Ti . 

In order to supplement the aforementioned nonrelativistic data at higher densities, we have also performed calcu- 
lations using the Jancovici p^ ] zero-temperature dielectric function. The results are shown in Fig. ^ for > 1 and 
rs = 0.0625 {x = 0.224), 0.0256 {x = 0.545), and 0.008 {x = 1.75). A comparison for the case of r, = 0.0625 confirms 
that the zero-temperature approximation works well at small Ts and large F (where 9^1 due to the relation (^); 
this conclusion is corroborated by inspection of Fig. |[ 

Thus the numerical results cover all values of and F (i.e., p and T) that are relevant for liquid EIPs. At Ts > 1 
and Fi > 1, the formation of bound states sets in. At ^ 10"^ and F^ ^ 1, the temperature reaches the values 
T > 3xWZ^/^ K, where the electron screening effects are completely unimportant. Finally, at F^ ^ 170, solidification 
takes place. 



B. Analytic formulae 

The calculated values of the screening free energy are fitted by the following function of Vg, Fg, and Z: 

f cdhVT\ + CTFa^e9i{rs)hi{x) 

' l+[bVT: + ag2irg)r-/rg]h2ix)' ^ ' 



The parameter 

Z r 



(1 + Zf/"^ - 1 - Z^/^l (30) 



ensures transition to the DH value of the excess free energy of the EIP, — —Z [(1 + Z)/3]^/'^ Fe^^ at small Fg. 
The parameter 

CTF = Coo^'/' (l - Z-^/^ + 0.2 (31) 

determines the screening in the limit of large Fg and small r^. The parameter Cao = (18/175)(12/7r)^/'^ = 0.2513 is 
consistent with the Thomas- Fermi approximation which becomes exact at small and very large Z (cf. Ref. 
|p6|). The parameters 



a= l.llZ"-475, 

& = 0.2 + 0.078 (lnZ)2, 

V = 1.16 + 0.08 InZ 



provide a low-order approximation to fif, (with a maximum error up to 30% at large Z and rs ft 1); while the functions 



9i{rs) = 1 + 
32(r.) = 1 + 



0.78 



21 + T,{Z/rsf 
Z -I 
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1 



1/2 



0.001 Z2 + 2re/ l + 6r? 



(32) 
(33) 



improve the fit at relatively large and reduce the maximum fractional error in /^e to 4.3%, and the root-mean-square 
(rms) error to ~ 1.5%. 

The factors hi{x) = [1 + {vp / c)^ Z~^/^]~^ (where vp = cx/^/l + is the electron Fermi velocity) and /i2(x) = 
(1 + a;2)~^/2 are relativistic corrections and may be omitted at a; ^ 1. 

Note that fie constitutes only a part of the ion excess free energy fa + fi^. The fit to this latter quantity is given 
by the sum of Eqs. ( p8|) and (|9|). The second and third columns of Table ^ present the rms and maximum relative 
differences between the calculated and fitted values of fa + fie at each value of Z. The comparison has been done 
for the set of finite-temperature numerical results at 0.1 < < 170 and 0.0625 < rs < 2.074. The remaining four 
columns of the table present the rms and maximum relative differences for the (ii + ie) internal energy and pressure, 
derived from the fits by the use of the thermodynamic relations 



Of 
dlnV 



PP/n 



Of 
Sin?', 



(34) 



The bottom line of the table is given for reference and presents the difference between the fit (^8|) and numerical 
HNC data in the OCP model (i.e., without the ie contribution). 

The calculated and fitted values of fie are shown in Fig. ^ for Z 6 and in Fig. ^ for Z — 1,2, and 10. For 
comparison, we have plotted the fit of Yakovlev and Shalybkov |l^] (YS) to their relativistic calculations, carried out 
in the zero-temperature approximation (justified at small rg and large Fg). In Fig. ^we have also plotted fie given by 
an analytic expression of Ebeling et al. reproduced by SB For the hydrogen plasma {Z = 1) it reproduces 
the Fade approximations of Ref. |^^. One can see that the fit of YS, in the range of its validity, agrees with our 
results. On the contrary, the approximation of Refs. [^,^ is clearly invalid in most cases. It exhibits unphysical 
behavior around F^ ~ 0.1, predicting an enhancement of screening with decreasing rg (e. g., for Z = I and F^ = 0.1 it 
gives larger fie at — 0.41 than at ~ 1.464). Moreover, the extrapolation to Z > 1, proposed in Ref. severely 
underestimates the screening effects. 

TABLE II. Root-mean-square and maximum relative differences between the fit and the HNC calculations for fu + fu, 
Uii + tiie, and Pa + Pi^; bottom line correponds to the OCP model. 



iSf/f) (%) 
rms max 



{5u/u) (%) 
rms max 



{5P/P) (%) 
rms max 



1 


0.6 


1.9 


0.9 


1.8 


1.2 


4.5 


2 


0.4 


1.1 


0.7 


1.8 


0.7 


3.0 


3 


0.4 


0.8 


0.9 


2.3 


0.7 


1.8 


4 


0.5 


1.2 


1.3 


3.1 


0.9 


1.9 


5 


0.6 


1.5 


1.6 


3.9 


1.2 


2.4 


6 


0.7 


1.8 


1.8 


4.3 


1.5 


2.8 


7 


0.6 


1.7 


1.9 


4.6 


1.7 


3.5 


8 


0.7 


1.5 


1.9 


4.6 


1.9 


4.1 


10 


0.6 


1.2 


1.8 


3.9 


2.0 


4.1 


12 


0.5 


1.2 


1.6 


3.2 


1.9 


4.5 


14 


0.5 


1.2 


1.9 


4.5 


1.6 


3.2 


20 


0.5 


1.1 


1.2 


2.8 


1.8 


4.5 


26 


0.6 


1.7 


1.2 


2.6 


1.3 


3.8 


OCP 


0.6 


0.7 


0.6 


0.8 


0.6 


0.8 
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FIG. 5. Calculated (filled circles) and fitted (solid lines) excess free energy, /cx ~ Xi{fii + fie) +a;e/ee, for Z — 1 for different 
values of the density parameter Vs- The dot-dashed line is the DH formula while the dashed and dotted lines represent the 
approximations of SB and IIT respectively. 

Figure ^ exhibits an analogous comparison for the excess free energy (^0|), for Z = 1. We have also plotted the 
Pade-approximation of IIT . Although the disgressions between the fit and numerical data of Ref. ||ll | lie within 
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0.4%, there are significant deviations between the IIT fit and our present results. This discrepancy originates from the 
relatively small number of numerical calculations used by IIT (32 computed values at 0.1 < < 10 and 0.1 < 9 < 10). 
Our fit, based on a much larger set of numerical data, not only reproduces these data but also the numerical results 
of IIT ]ll|. 
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FIG. 6. Calculated (filled circles) and fitted (solid lines) excess pressure of ions in the compressible electron background. 
Pi = Pa + Pie, in units of PoTi, where Po = UikBT. For comparison, the SB approximation is shown by dashed lines, and the 
DH and OCP approximations by dot-dashed lines. 
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FIG. 7. Equations of state (EOS) of fully ionized plasmas of four elements (Z — 4, 8, 13, and 26) given by the present 
analytic approximations. Solid lines show the pressure P vs. density p with account of the nonideality effects; dots represent 
EOS of the perfect gas of ions and electrons; dashes display an EOS in which the electron-ion screening effects are neglected. 
The gaps in some isotherms indicate the regions where the formation of bound electron states can be expected. 
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The excess ionic pressure Pa + Pie is shown in Fig. Calculated data are compared with the pressure obtained by 
differentiation (^4|) of our fit and of SB fit. The DH approximation, shown for reference, is calculated as the difference 
between the DH pressures of the electron-ion TCP and the electron OCP. The importance of the screening effects is 
verified by comparison of our calculated and fitted pressure with the pressure of the OCP in the rigid background. 
Pa, also shown in the figure. 

Figure demonstrates the validity of the EOS derived from our analytic formulae. The EOS of the perfect ion- 
electron gas is compared with the EOS which includes the nonideality of the electron and ion fluids but neglects the ie 
interactions; solid lines show the complete EOS. The gaps in some isotherms indicate the regions where the formation 
of bound states could not be neglected. The significant deviations of the broken lines from the full lines in certain 
ranges of p and T demonstrate the importance of the ion-electron screening effects. 



VII. MULTI-IONIC MIXTURES 



The multi-ionic mixture is a straightforward generalization of the previous single-ion model. In that case the 
effective Hamiltonian (||) reads : 

1 47re2 



2V ^ 

k#0 



e{k) 



(35) 



where p^k = ^ipn^ are the Fourier components of the ion charge number fluctuations. 

For the binary ionic mixture in a rigid electron background {e{k) = IV/c), the excess (non-ideal) free energy of 
the mixture, as well as the related thermodynamic quantities, can be expressed with high accuracy by the so-called 
"linear mixing rule" (LMR) in terms of the free energy of the pure phases : 

/cx(2'i, Z2,re,a;i) w a;i/cx(ri,a;i = 1) + (1 - Xi)fc^iT2,Xi = 0), (36) 
where F^ = TeZf^^ and xi = Ni/{Ni + N2)- The very high level of accuracy of the LMR (^6|) was first demonstrated 



by Hansen et al. |2l[] and confirmed later on by several authors, using very accurate MC calculations (e.g., pj3l|). 
The validity of the LMR in the case of an ionic mixture immersed in a responsive flnite-temperature electron 



background, as described by the Hamiltonian (35), has been examined by Hansen et al. | |21| in the flrst-order ther- 
modynamic perturbation approximation, and more recently by Chabrier and Ashcroft [ ^2| , who have solved the HNC 
equations with the effective screened potentials for arbitrary charge ratios ranging from a symmetric case (ZijZx ~ 1) 
to a highly asymmetric case [Z^jZi ^ 1). These authors found that the LMR remains valid to a high degree of 
accuracy when the electron response is taken into account in the interionic potential, except possibly for highly asym- 
metric mixtures in the region of weak degeneracy of the electron gas (where the departure from linearity can reach a 
few percent). 



VIII. CONCLUSIONS 



We have developed a completely analytic model for the free energy of fully ionized electron- ion Coulomb plasmas. 
The ideal part of the free energy of electrons and ions is described by Eqs. (^-(|l^) and is accurately represented by 
the analytic fits given by Eqs. ([l^)-(|2j). Note that these formulae provide the thermodynamic quantities of a free 
electron gas for any degeneracy and relativity. For the excess free energy of the electron fiuid at finite temperature, 
we adopt the analytic approximation from Ref. | [Tl[ |. For the excess free energy of the classical ionic OCP, we provide 
a simple interpolation (|2^) which accurately reproduces the Monte Carlo results at F^ > 1 and the Debye-Hiickel-Abe 
limit for F^ ^ 1. Finally, we have taken into account the ion-electron interactions by solving the hypernetted-chain 
equations for a large set of the parameters F^ , , and Z and constructed an analytic fit given by Eq. ( |29| ) . Our analytic 
formulae reproduce ja -I- /^e with accuracy ~ 1 — 2%, and the derivatives of this function with respect to Ts and F^ give 
an excess internal energy and pressure with relative errors not larger than a few percent. This analytic approximation 
is significantly more accurate than previous approximations of the free energy of the electron- ion plasmas. 

As mentioned in the introduction, our calculations imply full ionization, i.e. point-like ions from which their bound 
electrons are stripped completely. This model is realistic in various conditions at high temperatures or densities 
encountered in modern laser experiments and in various astrophysical situations like, for example, stellar, brown 
dwarf and giant planet interiors, or the envelopes of neutron stars. In these situations, complete ionization can be 
safely assumed. Furthermore, the present model can be used as the basis of more elaborated equations of state aimed 
at describing the thermodynamic properties of partially ionized plasmas and ionization equilibrium. Work in this 
direction is in progress. 
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